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Neutrino emission rate from baryonic matter in neutron stars via weak neutral vector interaction 
is computed up to order 0(y%), where vf is the Fermi velocity in units of speed of light. The vector 
current polarization tensors are evaluated with full vertices which include resummed series in the 
particle-hole channel. The neutrino emissivity is enhanced compared to the 0(v%) order up of 10% 
for values vf < 0.4 characteristic to baryons in compact stars. 



I. INTRODUCTION 

The pair-breaking neutrino bremsstrahlung is a impor- 
tant process contributing to the neutrino cooling of a 
compact star. The rate of these processes was computed 
initially at one-loop Q]-[3| and, more recently, including 
vertex corrections [j43 • Within the Standard Model the 
neutrino emission occurs via vector and axial- vector cur- 
rent interactions. The one- loop calculations suggest that 
the neutrino emission via neutral vector currents is large 
compared to the emission via axial vector currents. How- 
ever, the vertex corrections substantially suppress the 
emission via vector currents 0-0 while they leave the 
axial vector emission unaffected. As a result, the neu- 
trino emission via vector currents turns out to be sub- 
dominant compared to the same emission mediated by 
the axial- vector current interactions. Vertex corrections 
can be important in the context of neutrino scattering 
in proto-neutron-star matter, if the superfluidity sets in 
before the matter becomes transparent to neutrinos [8[. 

In this work we compute higher order corrections to the 
leading order (vanishing) contribution to the vector cur- 
rent neutrino emission rate. We adopt the same approach 
as in Ref. [5[, where a non-zero temperature propagator 
formalism was used to compute the vertex corrections 
to the weak interaction rates in pair-correlated baryonic 
matter. The convenient techniques of Nambu-Gor'kov 
propagators were used for a re-summation of particle- 
hole ladder diagrams - an approach first developed by 
Abrikosov and Gor'kov in electrodynamics of supercon- 
ductors Q (see also Ref. [10]). In this theory the re- 
sponse of superconductors to external probes is expressed 
in the language of propagators at non-zero temperature 
and density with contact interactions that do not distin- 
guish among the particle-hole and particle-particle chan- 
nels. It is equivalent to the theories initially advanced by 
Bogolyubov, Anderson and others, which are based on 
the equations of motions for second-quantized operators. 
A more involved approach was developed subsequently on 
the basis of the ideas of Fermi - liq uid theory for supercon- 
ductors/superfluids by Refs. pjJGS. The latter method 
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implements wave-function-renormalization of quasiparti- 
cle spectrum and higher-order harmonics in the inter- 
action channels, and postulates particle-hole (ph) and 
particle-particle (pp) interactions with different strength 
and/or sign. Application of this method in the context of 
neutrino emission and collective excitations can be found 
in Refs. [H, [T^j. In the preceding paper [E| the driving 
terms in both channels were taken to be identical (i.e. 
v pp = v p h) and equal to the lowest order Landau pa- 
rameter in baryonic matter. This identification for v p h 
is consistent with the theory of normal Fermi-liquids and 
guarantees the correct limiting form of the response func- 
tion in the unpaired state. The same identification for Vpp 
presumes pairing interaction which includes polarization 
of the medium in a manner that has been frequently used 
in the computations of neutron and neutron-star mat- 
ter pj|. 

It is now well established that the pair-breaking pro- 
cess vanishes to the leading order (LO) in small momen- 
tum transfer for both neutral vector and axial- vector cur- 
rents p]-[7j . The purpose of this work is to revise and ex- 
tend the results of Ref. [5] concerning the neutrino emis- 
sion via neutral vector currents. First, we recompute and 
correct the next-to-leading (NLO) contribution. Second, 
we extend the calculation to the next-to-next-to-leading 
order (NNLO) contribution and access the convergence 
of the series. 

The main focus and motivation of this work is the com- 
putation of the neutrino emissivity of baryonic matter in 
compact stars [HI, [l6| . The low-energy neutral weak cur- 
rent interaction Lagrangian describing the interaction of 
neutrino field ip and baryonic current j M is given by 

c w = -^=j^r(i-^, (1) 

where Gf is the Fermi coupling constant. The current 
of baryons for each 5-baryon is 

i„ = ^B7,.(4 B) -4 B) 7 6 )^, (2) 

( B) 

where ipB are the quantum fields of the baryons and Cy 

( B) 

and c A are the vector and axial vector couplings, re- 
spectively. 

The rate at which neutrinos are radiated from matter 
(neutrino emissivity) is given by (for derivations using 
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equilibrium or transport techniques see Refs. 0.E3) 



lm[L^( qi )U^(q)}5^(q-^2 qi ), 



(3) 



where qi = (<Ui,qi), i = 1,2 are the neutrino momenta, 
g(uj) = [exp(cj/T) — is the Bose distribution function, 
TLpxfa) is the retarded polarization tensor of baryons, and 



L^(qi,q2) = 4 + - ( qi ■ q 2 )g^ 



-It ^ QlaQ2/3 



(4) 



is the leptonic trace. Here the emissivity is defined per 
neutrino flavor, i.e., the full rate of neutrino radiation 
through weak neutral currents is larger by a factor Nf - 
the number of neutrino flavors. (We will assume Nf = 3 
massless neutrino flavors.) 

The paper is organized as follows. In Sec. [XT] we set 
the stage and introduce the correlation functions needed 
for a description of the superfluid baryonic matter. In 
SeaHIUwe discuss the vertex functions and polarization 
tensors for neutral vector current interactions of baryons 
and neutrinos and their double expansion in small pa- 
rameters of the theory. Section HVl computes the neutrino 
emissivity via neutral vector currents up to NNLO order. 
Section IVl contains a brief discussion of the results. Some 
details of calculations are relegated to Appendix [A] and 
the differences to Ref. [5| are discussed in Appendix iBl 



where F+ af (iuj ni p) = F aa , (iu n ,p), uj n = (2n + 1)ttT 
is the fermionic Matsubara frequency, cr y is the y 
component of the Pauli- matrix in spin space, u 2 = 
(1/2) (1 + (, p /e p ) and v 2 = 1 — u 2 are the Bogolyubov 
amplitudes, £ p is the single particle energy in the un- 
paired state corresponding to the momentum p, and 



e p = yj Q + A 2 is the single particle energy in the paired 

state, with A p being the (generally momentum- and 
energy-dependent) gap in the quasiparticle spectrum. 
The propagator for the holes is defined as G~^ a ,(iuj n ,p) = 
(5 crcr /(— iuo ni — p). For an £-wave condensate the spin 
structure of the propagators can be made explicit by 
writing G aa >(iuj n ,p) = 5 aa >G(iuj n ,p), F aa >(iuj n ,p) = 
—ia y F(iu n ,p), etc. 

The single-particle energies in the superconducting 
state, e p , are the solution of the Dyson- Schwinger equa- 

tion G~]r> = Ga,l',o ~ where G<t,<t',q is the diagonal 

free-particle propagator, whereas the self-energy Y, aj(7 t is 
a 2 x 2 matrix analogous to Eq. (J6j). The diagonal ele- 
ments of the self-energy matrix renormalize the quasipar- 
ticle spectrum and the density of states on the Fermi sur- 
face in the normal state. We will assume that such a pro- 
gram has been carried out and write the particle energy 
in the unpaired state as £ p = vp(p — Pf), where vp and 
Pf are the (effective) Fermi velocity and momentum. For 
contact pairing interactions ([5]) the gap function (the off- 
diagonal self-energy) is momentum-frequency indepen- 
dent; hereafter we set A„ = A. 



II. PROPAGATORS 



Consider low-density baryonic matter with attractive 
interaction in the ^o-channel. The interaction La- 
grangian for the baryons is 

tint = -V pp ^ ^B^(Pl)^B,i(P2)^B,l(p2)^B,t(Pl) 
Pi 

X S(p! +p 2 -Pi -P2), (5) 

where v pp < is the four-point coupling responsible for 
binding the particles in Cooper pairs. The imaginary- 
time momentum- space correlators are given by the 2x2 
Nambu-Gor'kov matrix 

V F^(tuj nj p) Gl a ,(iuj nj p) J 

The elements of the matrix are time-order correlators of 
the baryon field ipB and in the frequency- momentum 
domain these are given by [Ref. ||, Eqs. (5) and (6)] 

( u 2 v 2 \ 
G aa >(iuj n ,p) = 5 aa > h , (7) 

y ZUJ n - € p ZUJ n + € p J 

F aa >(iuj n ,p) = -i(JyU p v p (- — — ) (8) 



III. EFFECTIVE VERTICES AND 
POLARIZATION TENSORS FOR WEAK 
INTERACTIONS 



A. Vertex functions 



This section reviews and complements the discussion of 
the renormalization of the weak vertices in Ref. [5] . The 
diagrammatic approach and the mathematical structure 
of the theory is essentially the one found in Refs. 
UJj, but its physical content (the electroweak dynamics 
of baryons) is distinct. We start by stating the equations 
that are obeyed by the weak vector-current vertex = 
(To,r) in pair-correlated matter. We use the convention 
in which greek indices run through temporal and spatial 
1,2,3 components, i.e., \i = 0, 1, 2, 3. The latin indices run 
through spatial components only. 

The temporal (/i = 0) component of the vector-current 
vertex has four components in the Nambu-Gor'kov space 
in general. These components are given by the following 
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FIG. 1: A diagrammatic representation of the coupled integral equations for the effective weak vertices in superfluid baryonic 
matter. The "normal" propagators for particles (holes) are shown by single- arrowed lines directed from left to right (right to 
left). The double arrowed lines correspond to the "anomalous" propagators F (two incoming arrows) and F + (two outgoing 
arrows). The "normal" vertices Y\ and are shown by full and empty triangles. The "anomalous" vertices Y2 and Y3 are 
shown by hatched and shaded triangles. The horizontal wavy lines represent the low-energy propagator of Z° gauge boson. 
The vertical dashed lines stand for the particle-particle interaction v pp ; wavy lines represent particle- hole interaction v p h- 



integral equations [Ref. Eqs. (20)-(22)] 

fi = r + v p h(GTiG + FY 3 G + GT 2 F + FT 4 F), (9) 
T 2 = ^ pp (Gf 2 G + +Fr4G + + Gr 1 F + Ff 3 F),(10) 
f 3 = v pp (G + fsG + FTiG + G + T 4 F + Ff 2 F),(ll) 
f 4 = Y + v ph (G + Y±G + +FY 1 F + FY 2 G+ + G+Y Z F), 

(12) 

which are written in the operator form and are shown 
diagrammatically in Fig. [TJ The fi = index on vertices 
is suppressed, the bare vertex is To = 1, and the hats 
over the vertices indicate that these are 2x2 matrices in 
the spin-space, specifically 

The vertices Ti . . . T 4 are functions of the four- 
momentum transfer q and are independent of the incom- 
ing and outgoing momenta. Therefore, each term on the 
right hand side of Eqs. (f9l)- (fT2|) contains a polarization 
insertion defined as [Ref. [5], Eq. (25)] 



-i(j y TW and f( 3 ) 



n 



XX 



(<?) 



V M f 1 dx(X*X% 
4 J -i 



(13) 



where T is the temperature, X, X' G {G, F, G + ,F + }, 
p = (ipn,p) with p n being the fermionic Matsubara fre- 
quency, v(0) = m*pF/K 2 is the density of states at the 
Fermi surface, and x is the cosine of the angle formed by 
the vectors q and p. The convolution (or loop) is defined 
in Eq. ^ as 



+ q). (14) 



/oo 



Taking the integration limits over the infinite range 
— oo < £ p < oo is a valid approximation in the weakly 
coupled superconductivity. One may now exploit rela- 
tions among the loops [Ref. [5], Equations (26) and 
(A8)]: (G * F) ~ —(F * G), (G+ * G) ~ (G * G+), and 
(F * G + ) ~ -(G + * F) to show that Eqs. CEj) and CP 
have a solution only if T 2 +T 3 = 0, i. e., T 2 = — T 3 . Eqs. 
(j9j) and (fl2j) transform into each other on reversal of the 
time direction, which implies that Ti(a;, q) = 7T 4 (cg, g), 
where T is an operator affecting the reversal. It is +1 
for the scalar vertex and —1 for any vector vertex. This 
property is most easily checked by a calculation of the 
loops that change their sign, e. g. [F * F), for scalar 
and vector vertices. The solution of the linear system 
of equations for the temporal part of the vector-current 
vertex ©-CO]) reads 



r 2 = 



r c 



c 



v ph (A+C-BV+y 



c 



(15) 



(16) 



v ph (A+C 

where A r = U GG - n FF T, B = 2Ii FG , C 
(n GG+ + U FF )} and V T = U FG+ T+U GF . For T = 1 
one finds B = —T> + . Upon setting To = 1 we verify that 
the solutions (fT5j) and (fTB|l coincide with those given in 
Ref. [5], Eqs. (33) and (34). The vector vertex f is 
defined in full analogy to the scalar vertex: 



f 4 



r + VphiGTiG - FT 3 G - GT 2 F - Fr 4 F), (17) 
3 (Gf 2 G+ + FT A G + + Gf 1 F- FT 3 F), (18) 
,(G+f 3 G + FfiG + G+f 4 F - FT 2 F), (19) 

G+f 3 F), 
(20) 



= f + ^(G + f 4 G + -Ff 1 F-Ff 2 G+ 



(a) 



(b) 



(c) 



FIG. 2: The sum of polarization tensors contributing to the vector- current neutrino emission rate. Note that the diagrams 6, 
c, and d are specific to the superffuid systems and vanish in the unpaired state. 



where the same arguments as above imply f 2 = — I3 
and Ti = TT 4. There are two vectors available in the 
problem, therefore the most general decomposition of any 
vector vertex in terms of these vectors reads 



(21) 



where p = p/\p\, q = q/\q\ and Pi and Qi,2 are scalar 
functions, which depend on the moduli of these vectors. 
The solution of the set of linear equations (|17p - (fT2|) im- 
plies that 



Qi 
Q2 



= -Pi 



- v ph [A+C - BV+Y 
V- v ph lA + V - AD+] 



C - v p h[A + C - BT>+] 



(22) 
(23) 



where A = A~ (q • p) and V = V~ (q • p) are the first 
moments of the A~ and V~ integrals with respect to the 
cosine of the angle formed by the vectors q and p. 



B. Polarization tensors 



C. Expanding the polarization tensor 



We will use a double-expansion in small parameters 
which arise in the limit \q\ — >• with other parameters 
held fixed. Contrary to Ref. [5], we use the "symmet- 
ric" kinematics, where perturbed quasiparticle energies 
in the normal and superconducting state are written, re- 
spectively, as £± = £ p ± (v • q)/2 and e± = \J i± + A 2 . 
The loops ^4, S, C and V are products of the convolu- 
tion (F * F) and prefactors which depend only on the 
four-momentum transfer [Ref. [l2], Eqs. (15) and (17)]. 
These pre-factors are expanded in the small parameter 
S = \q\vF/oj. The expansion of the convolution (F * F) 
around £ p ~ £_ ~ £ + is carried out with respect to the 
second small paramater 77 = vp \q\/£,p] for details see Ap- 
pendix [XJ Specifically, we write the Taylor expansion as 



(F * F+) = ^(F * F+) n 5 n x n , 



(27) 



The 00 component of the vector polarization tensor is 
given by the sum of four terms 

n 00 = ToiGT^ + Ft 3 G + Gt 2 F + FT 4 F), (24) 

which are shown diagrammatically in Fig. [2j Exploiting 
the relations Ti = Tr 4 and T2 = — Ts one finds IIoo = 
r (^ + ri+Sr 2 ). Then, substituting the vertices ^ and 
([16]), we obtain 



n, 



00 



T (A + C-BV+) 
C-v ph (A+C-BV + ) 



Because B = — V + for T = 1 one recovers the polariza- 
tion tensor given by Eq. (35) of Ref. For the spatial 
components of the vector current polarization tensor one 
finds 



C - v ph [A+C - BT>+Y 



PiV F B 



V - v ph [A + V + AV+] 



Q31 



C-v ph [A+C-BV+) 
where we approximated the value of the P\ by vp- 



(26) 



where x = q • and the factor r] n /S n is absorbed in 
the definition of the coefficients (F * .F + ) n of the Tay- 
lor expansion ([2?]) . Substituting the expressions for the 
functions *4 + , 6, C, and P + in Eq. ([25]) we obtain the 
scalar polarization tensor truncated at order 0(rf) and 
0(J 6 ): 



(25) U 00 (q) 
2z/(0) 



xd 



1 — xS 



(F*F+) 



(F * F+) 



(l-x 2 £ 2 )(F*F+) 
4(5 4 



(l + x5)(F*F+) 



45 V 21 y v ;u 45 v ; 



1 



o(s 8 ), 



(28) 



where (...) = (1/2) Jj^ ..). Similar expansion for 
the spatial part of the vector current polarization tensor 



flU gives 

n»( g ) 

2v 2 f p(0) 



xS 




l-xS_ 


(x + x 2 5)(F *F+) - 



(l-x 2 5 2 )(F*F+) 

2S 2 22S 4 74S 6 

1 i- 

9 

U5 4 



135 
286(5 6 



x 2 5(F*F+) 
(F*F+) 



135 2835 



567 

(F*F+) 2 +0(£ 8 ). 



(29) 



In this last expression we need to keep only the terms up 
to £ 4 , because the vector vertices are proportional to the 
Fermi velocity, i. e., a small parameter. 



IV. VECTOR CURRENT NEUTRINO 
EMISSIVITY 




To obtain the neutrino pair bremsstrahlung emissiv- 
ity we contract the baryonic polarization tensor with the 
trace over the leptonic currents, see Eq. (|3j). The result 
can be cast as 



G 2 c 2 v N f 
48tt 4 



dujg(uj)ujJ(uj), 



(30) 



where cy = 1 for neutrons and cy = 0.08 for protons, 
Nf = 3 is the number of neutrino flavors in the Standard 
Model, and 



J(w) 



Jo 



dqq 2 (q 2 - o; 2 )Im [U 00 (uj, q) - IUi(w, q)] 



8u 5 v(0)v F 
405 



Tm(F*F+) [l + >yvjr], 



(31) 



where q is the momentum transfer and the coefficient 7 
is defined by Eq. (jA17|) of Appendix [XJ Substituting this 
result in Eq. ([30]) we find 



16G 



2^2 



e = 



12157T 3 



(32) 



where z = A/T and 

dy y 5 



I(z) 



vV 



1 



f(zyY 



41 

+ ) v 



(33) 



To order v F this result agrees with those obtained by 
Refs. 0, Ej. The numerical result at order 0{v F ) and 
their comparison to order 0(y F ) result are shown in 
Fig. [3j One observes that the corrections to the emis- 
sivity are small, i. e., the series expansion is justified. 
The emissivity turns out to be enhanced by about 10% 
for vp = 0.4 and by about 20% for vf = 0.6 in the tem- 
perature range < r < 0.92, where r = T/T c with T c 
being the critical temperature of pairing phase transition. 
In the close vicinity of T c the emissivity is reduced. 



FIG. 3: (Color online) Upper panel: Dependence of the in- 
tegral (|33|) on the reduced temperature r for vf = (solid 
line), 0.4 (dashed line) and 0.6 (dashed dotted line). Lower 
panel The ratio R of the integral (|33p computed at 0(v%) 
order to its value computed at order 0(v%). 



V. CONCLUDING REMARKS 

We obtained the neutrino emissivity via neutral vector 
currents in a perturbative expansion up to order 0(vp). 
The order 0{v%) term corrects our previous NLO re- 
sult [5[ and is in agreement with the expressions given 
in Refs. @, Q. (The expansion of the vector polarization 
tensor in Ref. [5] gave a mathematically spurious NLO 
contribution that is linear in the nucleon recoil [19|. As 
a result, the emissivity was suppressed by a factor T/m, 
which arises from the recoil term, instead of the correct 
factor Vp). 

From a practical point of applications in astrophysics, 
the 0(v F ) order term in the vector current neutrino emis- 
sivity is unimportant. Indeed, it is negligible compared 
to the axial vector neutrino emissivity as computed in 
Refs. p], [EEl, which remains unaffected by the vertex cor- 
rections. Nevertheless, the present reassessment removes 
the disagreement between Refs. 0, [H, and Ref. @, 
and confirms that the nonvanishing contribution to the 
vector-current emissivity arises at order 0(v F ). Simi- 
larly, it revises the results of Ref. [18| concerning the 
density response functions of cold neutron matter. How- 
ever, we note that the absence of 0(y F ) terms in the 
vector response is not protected by any symmetry of the 
theory (i.e., conservation law), whereas the 0(1) contri- 
bution remaining in one-loop calculations are prohibited 
by baryon number conservation. 
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The computation of the 0(v F ) order contribution to 
the vector-current neutrino emissivity shows that the cor- 
rections to the leading non-zero term are below 10% for 
values vp < 0.4 characteristic of baryons in compact 
stars. This result provides an evidence of the convergence 
of the series expansion of the vector-current polarization 
tensor in the regime where the momentum transfer is 
small compared to other relevant scales. 
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Appendix A: Expansion of the loop function 



The purpose of this appendix is to determine the coefficients of the expansion of the (F * F) loop in powers of the 
r] parameter, Eq. ([27]) . The Matsubara sum for the FF loop gives (Eq. (A12) of Ref. |5|) 




[/feO -/(**)] 



[l-f(e p )-f(e k )} 



(Al) 



For contact pairing interactions u p UkV p Vk = A 2 /4e p e/c. After an analytical continuation with retarded boundary 
condition, iu n — » uo + i5, we obtain for the convolution derived from Eq. (|A1|) 



(F*F+) = -— d£ p 



(fg ~ ffe) 



f(e p )-f(e k ) 



u 2 - (e p - e k ) 2 + iS 



(fg + ffc) 



1 - f(e P ) - f(e k ) 



w 2 - (e p + e fc ) 2 + iS 



(A2) 



To obtain the series expansion approximation to the inte- 
gral on the right hand side we expand first the integrand 
in series and carry out the integrations term by term. In 
doing so we will need to compute the imaginary parts of 
integrals of the type 



In 



f 

J — ( 



fn(e) 



(w 2 -4e 2 )" + iS' 



(A3) 



where at the order we are working we need only inte- 
grals with n = 1,2. For n = 1 the imaginary part of 
the integral is obtained with the use of the Dirac iden- 
tity. Upon the following change of variables £ — >> e = 
+ A 2 , £d£ = ede, A < e < oo, the integral ac- 
quires a factor of 2. A straightforward computation gives 

AGO 



Im/i 



Im 



7T 0(lj - 2A) 



2 v^r 

Consider now the case n = 



4e 2 + iS 



Im/o = Im 



y cx) 



4A 2 ' 

: 2, i. e., the integral 

/2(e) 

(o;2_4 c 2)2 +f< J 



We first change the variable £ — )• e 
ede and rewrite the integral as 



Im/o 



21m 



ede 



/ 2 (e) 



"A 2 (w 2 - 4e 2 ) 2 + id 



(A4) 
(A5) 

(A6) 



It is convenient to carry out a second transformation of 
variables defined as z = 4e 2 , dz = 8ede, 4A 2 < z < 
oo. Implementing the transformation we obtain 



ImJ 2 



2tt 



dz : 

4A 2 4y z - 4A 2 



f 2 (^/2)S { \z-oj 2 ) 



where we used the formula 
1 P 



D k+1 + id D k+1 



k\ 



(A7) 



(A8) 



where P denotes the principal value and 5^ (D) denotes 
the k-th derivative of the delta function. In the case 
k = 1 we obtain 



D 2 + iS D 2 
We further use the formula 



= ^+in5^(D). 



f(x)8'{x — a) = 
to obtain the final result 



Imio 



Im 



di 



>-2A)-| 



/ (x)S(x — a 

(oj 2 - 4e 2 ) 2 + id 
- 4A 2 



(A9) 



(A10) 



.(All) 
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To carry out the expansion we write the variables as 

£±=£(1±^), 77=^r, (A12) 



and substitute (]A12|) into the spectra e± = y £± + A 2 . 
Expanding the kernel to leading order we obtain 



)M4) 



Using Eq. ()A4|) one finds 

Im F*F+ o = / 2 ,A 2 tanh (^F 

^ V^ 2 — 4 A 2 



where e = y^ 2 + A 2 . The coefficient of the second order 
term is given by 



Im(F*F + ) 2 = 



ttA 2 0(uj - 2A) 



[hi (cj) + h 2 H], 



where the contributions arising from integrals wit n = 1 
and 2, respectively, are 



hi 



1 



rsech' 



U T) 



2T 



{2y 2 - 3) sinh (%=) 



4T 2 (2/ 2 "l)tanh (JL) 



1 27V 
(A15) 



h 2 = T-rrsech 2 f — ^ 



4^2 



V4T/ 



^- -Mzfi sinh (ii) 
2T y 2 -l \2TJ 



(A16) 



Note that the function hi )2 contain a factor rf/5 2 = 
(4y 2 )/(y 2 — 1), with y = A/2uj. The ratio appearing 
in Eq. ([33]) is then given by 



Im(F*F+) 2 1 



Im(F*F+) 4t/ 2 v4T/ 



h (If) csch (~~^ 



4T/ 



T y 2 



3) sinh 



f-1 

\2T J 



UJ 



(A17) 



Note that the second term in the brackets changes the 
singular asymptotics of the integrand in Eq. ([33]) at the 
end point y = 1 from (?/ — l) -1 / 2 to (2/ — l) -3 / 2 . 



Appendix B: Comparison with Ref. [Bj] 

In this appendix the equations from Ref. [5j are pre- 
ceded by roman I. As already mentioned above the po- 
larization tensor ([25]) coincides with Eq. (1.35). We now 
proceed to examine the three loops appearing in these 
equations. Substituting Eqs. (I.A1) and (I. A3) into the 
expression for the A(q) loop and using the definitions of 
the Bogolyubov amplitudes appearing after Eq. ([7]) we 
obtain 

A[q) = Swf [(e + e ' )(ee ' " ^' + A2) + w(? ' e " ?e ' )] * 

(Bl) 

with the short-hand notations £ = £ p , £' = £ p + q (asym- 
metrical kinematics) or £ = £ p _ q / 2 , £' = £ P +q/2 (sym- 
metrical kinematics), e = a/£ 2 + A 2 , e' = \J (£ 7 ) 2 + A 2 
and H3 



/(*) - /(O 



2ee' [ uj 2 - (e + e') 2 



(B2) 



On the other hand, substituting Eqs. (1.31) and (1.32) in 
Eq. (1.27), we see that the resulting expression does not 
contain the term uj{fc'e — £e') in Eq. (|B1|) , which vanishes 
in the limit gr — ^ 0. Using Eq. (I.A2) we further obtain 



B(q) 



2A J^[u;e' + (e' + e)eW. (B3) 



Substituting Eq. (1.31) in Eq. (1.28) we see that the term 
(e' + e)^' in Eq. (]B3|) is missing. Finally, using (I. A3) and 
(I.A4) in the expression for the C loop we find 



d 3 p 



2/(e) 



(2tt) 3 1 2e 



(e + e ')(«' + A 2 ) 
(B4) 



and we see that the term — uo{^e' + £'e) is missing in 
Eq. (1.29). Both terms that were dropped vanish in the 
limit qr — >• 0, because they are odd in £ while their convo- 
lutions involve integrations over symmetrical in £ limits 
[see Eqs. ([T3]) and Eqs. ([14]) ]. Thus, we conclude that, 
while the q — >> limit of the ^4, S, and C loops were cor- 
rectly identified, their small q expansion was carried out 
in Ref. [5| starting from incomplete expressions. 
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